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Abstract 

We compute on the lattice the "bag" parameters of the five AB = 2 operators of 
the supersymmetric basis, by combining their values determined in full QCD and in 
the static limit of HQET. The extrapolation of the QCD results from the accessible 
heavy- light meson masses to the i?-meson mass is constrained by the static result. The 
matching of the corresponding results in HQET and in QCD is for the first time made 
at NLO accuracy in the MS(NDR) renormalization scheme. All results are obtained 
in the quenched approximation. 



PACS: 11. 15. Ha (Lattice gauge theory), 12.38.Gc (Lattice QCD calculations), 13.75Lb (Meson-meson inter- 
actions), 14.40. Nd (Bottom mesons), 12.39. Hg (Heavy quark effective theory) 



1 Introduction 



This paper is devoted to a combined analysis of the matrix elements of the complete set of 
AB = 2 operators, which we computed on the lattice in both the static limit of the heavy 
quark effective theory (HQET) and in standard lattice QCD (with Wilson fermions). All 
five operators enter the phenomenological analyses of supersymmetric (SUSY) effects that 
might affect the Standard Model (SM) expectations for Am^ and/or Am^. It is therefore 
convenient to work in the so-called SUSY basis of operators: 



Ox = 




-7 5 )g*6- ? 7 / ,(l -7 5 )g J 


o 2 = 


¥{l - 


75)^^(1-75)^', 


3 = 


l\l - 


75)^(1-75)^ , 


4 = 


b\\ - 


75)^(1 + 75)^', 


5 = 


V{\ - 


75)5^^(1 + 75)^, 



where the superscripts denote colour indices, and q stands for either d- or s- light quark 
flavour. The first of the above operators has been widely studied over the last decade, since 
it is crucial for the SM description of the B° — B° mixing amplitude, whereas O2 and O3 
were also recently studied because they are relevant for the SM estimates of the relative 
width difference in the neutral 5-meson system, (AT/T) B 0. 

It is customary to parameterize the matrix elements of the operators ([!]) in terms of the 
so-called "bag" -parameters, which are introduced as a measure of the mismatch between 
the vacuum saturation approximation (VSA) and the actual value for each of the matrix 
elements, namely [fj, |3| 

(B a |di(/i)|B,°) = l^ljlBM, 



The hat symbol denotes operators renormalized in some renormalization scheme at the 
renormalization scale fi. To determine the values of the "bag" parameters £?x_s(^), we have 
performed a numerical simulation of (quenched) QCD on the lattice. While such a simulation 
can be made directly for the c-quark mass or somewhat heavier, present limitations of 
computational resources do not allow for a direct study of the 6-quark. For this reason 
we work in the range of heavy-light pseudoscalar masses mp G (1.7,2.4) GeV, from which 
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we have to extrapolate to the physical point, rriB d = 5.28 GeV and/or = 5.37 GeV. 
Although guided by the HQET scaling laws, this extrapolation is the dominant source of 
the systematic uncertainty in final results. To get around this problem we also computed 
the same matrix elements in the static limit of HQET on the lattice, and used them to 
constrain the extrapolations towards the physical point, rriB s/d - It is technically challenging 
to combine results from two different theories, namely one should match the S-parameters 
obtained in QCD onto the HQET ones so that the heavy quark scaling laws can be safely 
used. A special care to this issue will be given in the body of this paper. 
The main features of this work are: 

- The 5-parameters which appear in eq. (0) are computed using lattice QCD with 
Wilson fermions and are renormalized non-perturbatively in the (Landau) RI/MOM 
scheme. It is important to stress that we incorporated the recent proposal to remove 
the effects of the Goldstone boson contamination [|], |5| ; 

- The I?-parameters computed in the static limit of HQET on the lattice are matched 
onto the continuum MS(NDR) renormalization scheme. This matching has been made 
by using the one-loop (boosted) perturbative expressions || [7| ; 

- The conversion of the operators computed in lattice QCD from RI/MOM to MS(NDR) 
scheme is made at NLO accuracy [§, [|. Matching of the QCD operators onto the 
HQET ones has also been performed at NLO accuracy in a specified MS(NDR) scheme. 
In that procedure we use the recently computed 2-loop anomalous dimension matrices 
in HQET [Hj. With this matching at hand, we were able to constrain the extrapola- 



tion, i. e. to interpolate to the physical 6-quark mass; 

- Final results are presented in the RI/MOM scheme and in the MS(NDR) scheme of 
ref. ||. In addition, the parameters -Bi,2,3 are also given in the MS(NDR) renormal- 
ization scheme of ref. [|ll|]. [] 

The complete list of results can be found in table [TJ. Notice that we do not observe any 
SU(3) breaking effect in the S-parameters, i.e.: 



B 



B if) 



|o.99(2), 1.01(2), 1.01(3), 1.01(2), 1.01(3)} . (3) 



This paper is organized as follows: In sec. |2| we give the essential details of our lattice 
calculations and present the results as obtained for each heavy quark that we were able to 
access from our lattice and also in the static limit of the HQET. In sec. [| we outline the 
general strategy to combine the results of the two theories. We then explicitly give all the 
necessary anomalous dimension matrices and present the results of the combined analysis 
for all the five 5-parameters. In sec. |] we discuss the systematic uncertainties which are 
included in the results given in table [I]. We briefly conclude in sec. |[ 

Preliminary results were presented in ref. ]T2[. 
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Scheme 


RI/MOM 


MS(NDR) § 


MS(NDR) p 


B w 


(m 6 ) 


0.87(4) (if) 


0.87(4) (i|) 


0.87(4) (i|) 






0.82(3)(4) 


0.79(2)(4) 


0.83(3) (4) 


B$> 


(m&) 


1.02(6) (9) 


0.92(6)(8) 


0.90(6)(8) 


Bf 




1.16(3) (i 5 7 ) 


1.15(3) (i 5 7 ) 


- 




(m&) 


1.91(4) (I? 2 ) 


1.72(4) (if) 


- 


b[ s) 




0.86(2) (If) 


0.87(2) (^) 


0.87(2) (it) 




(m&) 


0.83(2)(4) 


0.80(1)(4) 


0.84(2)(4) 




(m 6 ) 


1.03(4)(9) 


0.93(3)(8) 


0.91(3)(8) 


b { : ] 


(m b ) 


1.17(2) (3) 


1.16(2) (i 5 7 ) 






(m b ) 


1.94(3) (If) 


1.75(3) 





Table 1: The main results of this paper: B -parameters defined in eq. renormalized at 
fj, = mi, = 4.6 GeV and in three renormalization schemes: RI/MOM, MS of ref. j^/ and the MS of 



ref. [11]. The results are obtained in the quenched approximation. 

2 Direct Lattice results 
2.1 Computation in lattice QCD 

In this subsection we recall the main elements of our lattice simulation, the details of which 
can be found in refs. |I3| , [14]]. We work with a lattice of the size 24 3 x 48, at /3 = 6.2, and 



use the non-perturbatively improved Wilson action \TE\. Note however that the 4-fermion 
operators, which are the main target of the present work, are not improved. Our data-set 
consists of 200 independent gauge field configurations. We work with 3 values of the heavy 
and 3 values of the light quark masses, corresponding to the Wilson hopping parameters: 
n q G {0.1344,0,1349,0.1352}, and k q G {0.125,0,122,0.119}. The mass spectrum and the 



decay constants have already been discussed in our previous publications |TJ, [TJ], [RJ and 
we immediately turn to the computation of the I?-parameters. 

The starting point is to compute the 2- and 3-point correlation functions 



c?kt) = {Y i j(x,t)ji<p)) 



X 



t>0 



Z, 



2 sinh Mj 



-Mjt 



3 



x,y 

T-t^o __y/Zp_ M ptl . (p\0-(u)\P) ■ V ^ e~ Mpt2 ( I ) 
2sinhM P WU t w\r q ) 2ginhMp e , 



where the operator Oj is placed at time equal to zero, the source of the pseudoscalar mesons 
(P5 = (375(7) is fixed at £1 = 16, while the other source operator moves around the periodic 
lattice (of the size T = 48). At some £2 = t, which is sufficiently far from the first source 
and from the operator Oj, the lowest lying heavy-light pseudoscalar meson, P g , is isolated. 
J in the above equations stands for either P 5 or A = Q7o75<?- 
To extract the parameter Bi(fi), one computes 

p U\ Cf 3) (tl,t;/^) T -t»0 (Pg\6l(fj)\Pg) _ , , , , 

-z\ 6%{t) cS(*i) -|<o|A m>l 2 

where the 2-points functions are used to eliminate the exponential terms from the 3-point 
functions and also to divide out the (8/3) fpm P from eq. (0), thus accessing directly the 
wanted P-parameter. Similarly, to reach other four P-parameters we form the ratios 

Rb ( t) = cf> ^ w im^> s b M , (6 ) 

6 4 Cf p (() C« (d) 6, \(0\PM\P,)V 

where 6j G { — 5/3,1/3,2,2/3} for z e {2,3,4,5}, respectively. 

As already discussed in our previous works [|T3j, |H|| , the operators Oi(n) are renormalized 
non-perturbatively in the (Landau) RI/MOM scheme by using the method explained in 
detail in ref. |T7| (see also references therein). The method is based on the possibility 



of computing amputated 4-quark vertices at sufficiently large quark virtualities with the 
operators inserted at zero momentum. The RI/MOM renormalization condition is merely 
imposed on various projections of the amputated Green functions, which then lead to a full 
set of 9 renormalization g^)) and 16 subtraction (A(#q)) constants, where g\ = 6/ ft is 
the bare lattice coupling. The renormalized operators are given by 

Oibi) = ZvQi,g*) [0 3 {gD + A jk (g 2 )O k (g 2 )]^. . (7) 

Since the computation of the off-shell quantities (4-quark vertices) is performed in the 
Landau gauge, the scheme is often referred to as the Landau RI/MOM scheme. 

In practice, one computes all the renormalization and subtraction constants with specific 
(nonzero) values of the light quark masses and then extrapolates all Z's and A's to the chiral 
limit. It has been pointed out in refs. |5| that such an extrapolation can be contaminated 
by the Goldstone boson (GB) contributions. The recipe to subtract these contributions 
away has been proposed and implemented in ref. |18|| (see also appendix of this paper). We 
employed that prescription and recomputed the renormalization and subtraction constants. 
Their values, together with Zp(p) and Za which are needed in eqs. (||,|), are listed in 
appendix. This is a new feature of our computation which improves (corrects) our previous 



results, presented in refs. |L3| . 
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Kq = 0.122; K q = 0.1349 

NPR (in RI/MOM) at 2.8(1) GeV 
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Figure 1: Signals for the ratios Rs^t), defined in eqs. (H) and are shown for the combination: 
K q = 0.1349, kq = 0.122. The operators are non-perturbatively renormalized (NPR) in the 
(Landau)RI/MOM scheme at u. ~ 1/a = 2.8(1) GeV. 



In fig. |T|, we illustrate the quality of the signals for the ratios i?s 1 _ 5 (t) for a specific 
combination of heavy and light quarks. After inspecting the ratios for all 9 pairs of kq-k q , 
we find that common stability plateaus are reached for 



Rbi 

Rb 2 , 3 



t E [28, 35] ; 
t E [30, 35] ; 
t E [29, 35] . 



(8) 



On each of these plateaus we fit the ratios to a constant and hence extract the parameters 
Bi^(ii) for each combination of the heavy and the light quark masses. Every parameter 
is then linearly interpolated in the light quark mass to the strange and extrapolated to the 
up /down light quark mass, by using the (standard) lattice-plane method (see refs. [0 for 
details). The extrapolations are very smooth as it can be seen from fig. |2|. In table 0, we 
present a detailed list of the results for all five 5-parameters renormalized non-perturbatively 
at three different values of the renormalization scale (i in the (Landau)RI/MOM scheme. 
The results are given for the values of the heavy quark masses, which are in the region of 
the charm quark and slightly higher, and for the light quark interpolated to the strange and 
to the averaged up/down mass. 
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Figure 2: Extrapolations to the light up/down and interpolations to the strange light quark 
mass are shown for all five B-parameters in the case of fixed heavy quark mass corresponding to 
kq = 0.122. Empty symbols denote the B-parameters directly accessed in our simulation. Filled 
symbols correspond to the B parameters extrapolated to the physical u/d and s quark respectively. 
In the figure [xa = 1.03. 



2.2 Computation in the static limit of HQET 

To avoid a confusion in notations, we consistently use "tilde" symbols over the operators 
and the .B-parameters computed in HQET. Instead of the five operators that we listed in 
eq. flU), in HQET one deals with only four of them, namely 

a, = /^(i- 75)^^(1-75)^', 

Oj = ^(1-75)^(1-75)^', 

4 = ^(1-75)^(1 + 75)^', 

5 = /i l (i- 75)^(1 + 7 5 )g\ (9) 
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q=up/down 


q=strange 


Scale (fi) 


1 QH "I HpV 


9 8f1 "l HpV 


O. C/ I *L J \J V V 


I.jI 1 1 VJC V 


9 8(1 "l HpV 


3 qf9^l HpV 




0.856(20) 


0.841(20) 


n 007 /on\ 

0.825(20) 


r\ z' 1 /' 1 / 1 7 N 

0.866(15) 


n orn/i r\ 

0.850(15) 


n o 7 / 1 7 \ 

0.835(15) 




0.875(38) 


0.817(31) 


0.797(31) 


0.904(27) 


r\ a / oo\ 

0.843(22) 


0.822(22) 




1 on c / to\ 

1.295(7o) 


1 nTo / c 1 \ 

1.072(51) 


n n /i 1 / .4 o\ 
(J. 941(42) 


1 OO /I / .4 O \ 

1.284(48) 


1 C\ r 7 A f a\ 

1.(J74(g4) 


0.947(28) 


B[ Ql \ri 


l.U92(2o) 


1.10/(2z) 


1 no c / on\ 
l.UoO(2U) 


1.112(lo) 


1 1 on / 1 ^ \ 
1.129(lo ) 


1.056(15 J 


B^V) 


1.523(37) 


1.322(29) 


1.332(27) 


1 £?n^ / on\ 

1.606(30) 


1.386(24) 


1.392(22) 


Bj Q2) (/i) 


r\ o o n /o 7 \ 

0.880(25) 


r\ o/" 1 O /O \ 

0.862(24) 


0.849(24) 


n 000/1 t\ 

0.883(17) 


0.866(17) 


n 7 / "i t\ 

0.852(17) 


Bf 2 V) 


0.886(37) 


f\ OO /I /OI \ 

0.824(31) 


r\ on o / n\ 

0.803(30) 


n n a n / 0/" \ 

0.940(26) 


n ^7 a ( \ 

0.874(22) 


n 071/ nn\ 

0.851(22) 


B^V) 


1.275(81) 


l.l)b(J(55) 


0.932(45) 


1 On /I / /I O \ 

1.294(48) 


1 non /oo\ 

l.(J89(oo) 


n nro / 00 \ 

(J.9oo(28) 


Bf 2 V) 


l.l)bo(2b) 


1 no /o t:\ 
l.Uoo(2o) 


1.01o(24J 


1.124(1 / ) 


1 1 A A ( 1 T\ 

1.144(1 / ) 


1 n 7n / 1 ^ \ 
l.U f U(lb ) 


Bf 2) (/i) 


1.673(45) 


1.427(36) 


1.421(34) 


1 TO /I / O £ \ 

1.784(35) 


1 C on /oo\ 

1.520(28) 


1 7 11/ O/ 1 \ 

1.511(26) 


B^V) 


0.887(26) 


0.869(25) 


0.856(25) 


0.886(16) 


0.867(15) 


0.855(15) 


Bf 3) (/i) 


O.y50(oo ) 


0.881(2 ( ) 


0.85< (26J 


0.y62(ZO J 


o.8yo(iy J 


o.86y(iy J 




1.349(84) 


1.123(54) 


0.989(43) 


1.288(51) 


1.090(34) 


0.966(27) 


Bf 3) (/i) 


1.114(23) 


1.136(21) 


1.062(20) 


1.119(14) 


1.141(13) 


1.068(12) 




1.885(47) 


1.592(35) 


1.573(33) 


1.925(32) 


1.622(25) 


1.601(23) 



Table 2: B -parameters, as defined in eq. (^), extracted from our non-perturbatively renormalized 
data at three values of the renormalization scale [i in the RI/MOM renormalization scheme. From 
top to bottom of the table, the bag parameters correspond to the heavy quark with kq x = 0.125, 
kq 2 = 0.122 and kq 3 = 0.119, respectively (they are separated by the horizontal lines). 



where h stands for the infinitely heavy (static) quark. In the HQET, the operator O3 is 
related to 0\ and 2 by the equations of motion as 



-o 2 - l6x . 



(10) 



The computation of the first two operators has been explained in detail in refs. || [2(J . The 
data-set consists of 600 configurations gathered on the 24 3 x 40 lattice at (3 = 6.0. The light 
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quarks were simulated by using the tree-level improved Wilson action and three values of 
«, E {0.1425, 1432, 1440}. 

To reach the HQET parameters equivalent to the ones appearing in eq. (0), which we 
will call Bi, one computes the following 2- and 3-point correlation functions: 



-A£t 



Cf A (t) = (J> (x,t)4(0)} ^ Z A e 

X 



x,y 



t>0 j? (-Pg|Q»MI-Pg) -A£(ti-t 2 ) Qix 

A ' 2M P ' 1 ' 

where the source operators are the axial currents which are extended by using the so-called 
cubic smearing procedure [21]. AS in eq. (|TT1), is the binding energy of the heavy meson P q . 



Note also that in this case, y Za = (0\A \P q ) / y/2Mp. To get the desired 5-parameters, we 
form the following ratios: 



n sM = g|^4 m «%M> s g iM , (12) 



where bi e {8/3, —5/3, 2, 2/3} for z <G {1, 2, 4, 5}. In this case, the fixed time has been chosen 
to be t 2 = 35. In fig. [3], we show the quality of the signals for all four Ri(ti). On the plateaus 
we fit to a constant and thus obtain the values of the corresponding parameters Bi(fi). Our 
data is renormalized in the MS(NDR) scheme, after matching the lattice regularization 
scheme onto the MS(NDR) by using one-loop boosted perturbation theory, as explained 
in great detail in ref. 0]. The matching scale fi = q* is varied between 2/ a < q* < 
7T /a, and the results are run to /i — = 4.6 GeV. The spread of values is assigned 
to the systematic uncertainty. As in the previous subsection, all 4 operators are linearly 
extrapolated (interpolated) in the the light quark mass to up/down (strange). The resulting 
values are listed in table [| In that table errors are statistical only, obtained by using the 
standard jackknife procedure. 



3 Extrapolation to the B-mesons 

Armed with "raw" results obtained in full QCD (table |2|) and in HQET (table |3|), we now 
discuss the extrapolation of the QCD results to the physical 5-meson mass. The aim of this 
section is to provide a consistent way to constrain that extrapolation by the static HQET 
results in order to reduce the systematic uncertainties. 

A common wisdom is to follow the HQET scaling laws, according to which every B- 
parameter scales with the inverse heavy quark (meson) mass as a constant, and to ex- 
trapolate to the desired heavy meson mass. The (unknown) 1/mp corrections are to be 
determined from the fit with our data. To use these scaling laws, however, one first need to 

2 In the static limit the axial current is identical to the pseudoscalar density. 
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t (Time of the moving source) t 



Figure 3: Signals for the ratios i^g^t) defined in eqs. from which the Bi parameters are 
extracted. The common plateaus are chosen for t\ = t G [32, 36]. The plotted ratios correspond to 
K q = 0.1432. The operators are renormalized perturbatively and evolved to /u = mj,. 



Light quark 


q=up/down 


q=strange 


B l (m b ) 


0.89(4) 


0.89(4) 


B 2 {m b ) 


0.82(4) 


0.83(3) 


B±(m b ) 


1.07(4) 


1.07(4) 


B 5 (m b ) 


2.37(10) 


2.40(10) 



Table 3: B -parameters, extracted from our HQET data by using the boosted perturbative 1-loop 
matching onto the continuum MS(NDR) renormalization scheme. 

relate the matrix elements of the QCD operators, (Oi(ri)) from eq. ([!]), to the HQET ones, 
(Oj(/i)) of eq. @ pi}] . This matching is made in perturbation theory at some suitably cho- 
sen renormalization scale, for example /i = m b . The matching is crucial since the anomalous 
dimensions for these operators in the two theories (QCD and HQET) differ. Moreover, when 
dealing with the 4-fermion operators, it is important that the matching between the two 
theories is made at NLO accuracy because it is at this order that the scheme can be fully 
specified (leading order anomalous dimensions are universal). Before entering the details of 



9 



that matching, we now outline the basic strategy that must be followed. 

Matching of the QCD operators, renormalized at some high scale > mp, and the 
HQET ones, renormalized at some low scale \il <C mp, is made at /i = mp by using the 
following expression Q 

W T QCD [m P ^}- 1 {6{^)) mp = C(m P ) W^^mp,//]- 1 (0(p'))+o(— ) + ... (13) 

where Wq C£) [/12, Pi\~ l is the matrix encoding the full QCD evolution from a scale \i\ to 
/i 2 of all five AB = 2 operators which are, for convenience, collected in a five-component 
vector (0(p)) mp . Likewise for W^g £T [/i 2 , A^i] -1 in HQET. These matrices will be speci- 
fied later on. We will be working in the MS(NDR) scheme in which the matrices of the 
anomalous dimension coefficients are known at NLO in both theories. Hence, the matching 
matrix, C(mp) = 1 + J2 n C ^[ a s( r mp)/An] n , is also known at NLO, i.e. is completely 
determined El. 



On the HQET side, we also consider five operators, (O (//)), where we add the operator 
O3 by means of eq. fllCf). In this way the matching matrix is squared (5x5). We now 
put all the evolution expressions appearing in eq. (0) on its l.h.s. 

w hqet[^', mp]' 1 C' 1 {mp)Wl CD [m P , ft]' 1 (6{fi)) mp = 

( ^N T Q cD{^P^]C{mp)W T HQET [m P ^r l Y \d(fi)) mp = (O(m0> + O (—) +... (14) 
\v v // \m P j 

so that the l.h.s. manifestly satisfies the HQET scaling laws which are the intrinsic prop- 
erty of the r.h.s. One proceeds similarly for the bilinear operators to define the matching 
constants M.^m P , fi, //]. In terms of £>-parameters, eq. fll~4] ) then reads 

Ml[m P , At, f/\ (h- l Mi l [m P , p, //]b) B(fi) = B(fi') + O (-^j + ... (15) 

where b is the diagonal matrix of the coefficients appearing in the definitions (0), i.e. 

b = diag(8/3, — 5/3, 1/3, 2, 2/3), and by B(fj) (£?(//)) we designated the vector column of 
our five 5-parameters. 

Based on the above discussion, a simple recipe can be applied to our data, namely evolve 
to the same \i = p! and create the quantity 

$(mp, fi) = M 2 2 [m P ,fi] (b^M^mp./ijbj • B(/i) (16) 

which can be fit either freely as 

$(m P , yU ) = a ( / i) + ^M, (17) 

mp 



3 Instead of matching at the point corresponding to the mass of the heavy quark, we choose to do it at 
the mass of the heavy-light meson, mp. 
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where ao(/x) and ai(/x) are the fit parameters, or by constraining it by the static HQET 
results, i.e. 

%^) = a» + *l + ^, (18) 

mp mp 

where the coefficient a' (fj,) is constrained by the static value, B(/j), so that one can probe 
the term 0(1/ mp). As a result of these two procedures, we obtain the HQET values of 
the 5-parameters, i.e. $(mg s/d ,/i), which are then to be matched back onto their QCD 
counterparts. 

To keep the expressions as short as possible, we will now split the discussion into two 
pieces: we will first discuss the extrapolation of the first three 5-parameters and then the 
last two. This can be done because all the matrices, Wq CD [/x 2 , ^hqetI^, A*i] and 
C(n) are the block- matrices of the form [3 x 3] © [2 x 2]. 



3.1 



Getting the physical results for 



,2,3 



At the leading order in perturbation theory, the anomalous dimensions for our 5-parameters 
in the RI/MOM and MS schemes are the same. This is not the case at NLO, and in 
order to proceed we need to convert our RI/MOM results from table ^ into the MS(NDR) 
scheme. It is crucial to specify the set of evanescent operators or the Dirac projectors used 
to renormalize the operators because only with this information at hand, the MS(NDR) 
scheme is unambiguously defined [p5|| . In this subsection, we will use the MS(NDR) scheme 
of ref. ]TTJ (see eqs. (13-15) of their paper) in which the Wilson coefficients for the SM 
expression for the (Ar/r)# s have been calculated at NLO. Therefore, the results for the 
5-parameters that will be presented in this subsection can be directly combined with the 
Wilson coefficients of ref. [OJ . 

The conversion of the operators Of!/^ 101 ^ (ji) and P 5 RI//MOM (/i) to the MS scheme is pro- 
vided by the following expressions 



(Ox(ri) 
{OM) 



MS 



I + ^123 



a s (n) 

47T 



<0i(A*)> 

<02(aO> 



RI/MOM 



{PM) 



MS 



1 + (P 5 (^)) RI/MOM , 



where the NLO matching coefficients are given by 

-42 + 72 log 2 



rp 



16 
~3~ 



^123 



1 

9 








m 



(19) 




61 + 441og2 -7 + 281og2 ] .(20) 
-25 + 28 log 2 -29 + 44 log 2 



To get the central values, we will convert the results from table |2|, obtained at fi — 2.8(1) GeV 
and run them to fi = m\, = 4.6 GeV [22]. It should be noted that the matching RI/MOM — > 
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MS(NDR), made in ref. [n|, was incorrect because the MS scheme was not the one of ref. [|TT|]. 
but rather the one of ref. 0. Although the numerical differences are very small, the physical 
results presented in ref. are not fully consistent because the matrix elements matched 
onto the MS scheme of ref. 101 were combined with the Wilson coefficients of ref. Ill 111 . The 



correct physical results were presented in ref. [23 



The evolution from the scale fi to mj, in this MS scheme, is described by Jl I 



Oi{m b ) 
2 {m b ) 
3 (m b ) 



MS 




MS 



(21) 



where the operator W QCD[ m bi f4 = M(m b )U(m b , /z)M -1 (yu) contains the information on the 
evolution obtained at the leading (t/(/i, m h )) and the NLO (M(/z)) in perturbation theory. 
For our purpose, it is convenient to write the evolution matrix in the following form 



W 



QCD 



m b , /i] = w(m b )w (/i) 



(22) 



where 



and j3 



Ufa) = M(/i) a a (ji)^' 2 *> f (23) 
11 — 2np/3. The scheme independent, one-loop anomalous dimension matrix is 



To 



whereas the NLO contribution 



M(/i) 





-28/3 4/3 
16/3 32/3 



4tt 



(24) 



is encoded in the matrix 



/ 



J 



MS 
123 



485 
242" 






V 



o 

_4592 

wm 

1612 
3TJ25" 





19083 
^025" 

36233 



(25) 



/ 



We have set np = 0, since our lattice results are obtained in the quenched approximation. 
The evolution of the pseudoscalar density is given by 



(P 5 K)) MS = 

where 7p = — 8 and Jp s 



-7P/2A, r 



a s (m b ) J 

998/363, for n F 



1 + 



a s (p) - a s {m b ) m 
— J p 



47T 



MS 



(26) 



0. 



12 





q=up/dowri 


q=strange 




u. 1Z0 


n 1 99 

u. 1ZZ 


u. 1 iy 


n 1 ok 

U. 1ZO 


n 1 99 

u. ±ZZ 


u. 1 iy 


m P [GeV] 


1 7^0^ 
1. 1 o^y j 


Z.UZ^IU J 


z.zo^iz j 


l.oO^o j 


z. iz^y j 






0.891(21) 


0.913(26) 


0.920(27) 


0.900(16) 


0.916(18) 


0.919(16) 




0.745(29) 


0.771(29) 


0.824(26) 


0.788(21) 


0.817(21) 


0.833(18) 




0.896(41) 


0.889(43) 


0.943(41) 


0.902(27) 


0.918(26) 


0.921(26) 



Table 4: B-parameters, extracted from our lattice QCD data after conversion to the MS(NDR) 



scheme of ref. [11] at u = 2.8(l)GeV and running to fx = m b . 



With all of the above formulae at hand, we convert our .B-parameters from table. [2| to 
the MS scheme (at /1 = 2.8(1) GeV), run them to \i — m b , and list their values in table (|. 

The next step is to match the I?-parameters from table f| from QCD to the HQET 
where we can use the heavy quark scaling laws to extrapolate each .B-parameter to the 
physical point, i.e. to mp, d/s - The crucial ingredient that enters the matrix A1j 1 [mp,mb] = 

•M ^[mp, m b , m b ] of eq. QI^), and hence the quantity l(mp, m h ) in eq. flTED, is the matching 
matrix C(mp) which relates the QCD operators, computed in the MS(NDR) scheme of 
ref. |Tl|], onto the HQET ones computed in the MS(NDR) scheme of ref. || (and vice 
versa). In ref. it has been shown that (at NLO) this matrix has the following form: 



rr (i) cxsjmp) . 
123 4tt ' 



"123 




61/12 
-77/12 




(27) 



where we introduced the index "123" since we consider only these operators here. 

The last piece of information needed to construct Ai^ 1 [mp,m b } of eq. (0), is the evo- 
lution operator in the HQET at NLO. In a notation analogous to eq. fl2"2|) we have 



(ON) 



™ T hqet 



m b ,u] 1 (0{fj)) 



W HQET [m b ,u] = w(m b )w 1 (fi) 



(28) 



where 



w{n) = M(fi) a s (/i)-^ /2/3 ° 



(29) 



In this case 10 



3 
To= ~ I 2 1 , 
" 1 1 2 



13 



¥ , T MS 

123 "17" 



(30) 



MS 
123 



317 68tt 2 



1089 297 






2639 7T 2 
"8712 + 99 
3181 74vr 2 
+ ~297 
2vr 2 



\ 



4356 
2639 



1913 2tt 2 



4356 
3907 



V 



4356 99 



+ 



99 
74vr 2 



4356 297 



(31) 



/ 



where, as before, we have set rip = 0. 

To obtain the quantity $(mp,m&) of eq. Q16D , one also needs .Mj hPiH- At NLO 



this information can be extracted from ref. |Pq| . The matching of the axial current and of 
the pseudoscalar density is given by 



An 



i 



3 4vr 



1 + 



8 a 8 (jn) 
3 4tt 



A)(^) , 
2 (/i) . 



(32) 



while the expression for the running of the axial current (for rip = 0) is 



A Q (m b ) 



a 



-7 A /2/3o r 



a s (m b ) J 



439 



287r 2 \ a s (/x) - a s (m b ) 



1089 297 



47T 



M»), (33) 



where the leading order anomalous dimension is given by 7^4 = —4. 

By combining all of the above ingredients, we form the quantities $i,2,3(^p, wii>), use 
eq. ( jTTD , and extrapolate to rriB d/s - The result of that extrapolation is: 



1.038(91) 
1.055(64) 
1.117(94) 



^1,2,3(^,^6) 



1.009(36) 
1.045(24) 
1.070(48) 



(34) 



where we wrote separately the results of the extrapolation of our data with the light quark 
extrapolated to d (left), and those with the light quark interpolated to s (right). By including 
the .B-parameters from table [3|, we fit our data to eq. (|18D , from which we get 



$1,2,3 {mB d ,m b ) 



0.976(42) 
0.943(28) 
0.845(52) 



$1,2,3 {m Bs ,m b ) 



0.967(21) 
0.955(17) 
0.857(31) 



(35) 



These two extrapolations, for each component of the vector $i,2,3(^p, Tn b ), are illustrated in 
fig. § for the case of the light <i-quark. Results of that extrapolation are the S-parameters 
in HQET, which we then match back onto their QCD values to get our final results for 
-B-parameters. These numbers are presented in table |5|. 
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1.2 
1.0 
0.8 
0.6 

1.2 
1.0 
0.8 
0.6 

1.2 
1.0 
0.8 
0.6 



— c|— 



•B,'(m b ) 



-f 



oijm^Jii] 
• B 3 (m b ) 



0.0 



0.5 



1.0 



1.5 



1/Mr 



Figure 4: Extrapolation to the physical Bd meson mass (squared symbols) in the inverse heavy 
meson mass. The unconstrained linear extrapolation for each of the components of the vector 
^i,2,^{fnp-,fnb) from our data (empty circles) to ^1,2,3 (^_B d , mb) (empty square) is depicted by 
the dotted line. The result of the constrained extrapolation (filled squares) by the static HQET 
bag-parameters (Riled circles) is marked by the dashed line. i th component of the vector is marked 
by[i\. 



An important issue to be mentioned is the treatment of the statistical errors when 
constraining by the static HQET results because the three points obtained in QCD are cor- 
related among themselves and are uncorrelated from the one obtained in HQET. Although 
it may look trivial, we prefer to mention how these errors have been treated in this work. Q 
For each jack of our QCD data we form the so-called augmented xl by including the HQET 

4 For a clear discussion about the treatment of the statistical errors in such a situation, see ref. p7|. 
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value Bj and its error aj as: 



a: 



0)3 



(Ah 

m P . 



m 



2)3 



Pi 



o-,-[$(m Pj )] 



+ 



<74 



(36) 



where "[7]" denotes the component of the vector $(mpj, the error of which is a~j[$(mp.)] 
for the i th value of our three heavy-light meson masses, mp i . By minimizing (Xa)j' we find 
the values of the parameters (a' 1 2 )j, where (a' Q )j is constrained by the prior knowledge of 
the static result. In this way we get the value for )[?'], for each jack of our QCD 

data. The final error, <Tj($(m Bd/ J), that we quoted in eq. fl35|), is obtained as a standard 
error over all jacks (Njk) 



\ 



Nj K -I 
N 



JK 



Njk 
k=l 



d/i 



N 



JK 



,fe=i 



• (37) 





unconstraiiied 


constrained 


Light quark 


d 


s 


d 


s 


5p(m fe ) 


0.938(81) 


0.905(32) 


0.875(37) 


0.867(18) 


5f(m 6 ) 


0.923(32) 


0.915(30) 


0.826(25) 


0.836(15) 


5p(m fe ) 


1.192(101) 


1.141(52) 


0.901(56) 


0.914(33) 



Table 5: Final results for the first three B-parameters defined in eq. (|), in the MS (NDR) scheme 
of ref. at \i = mj,. 

From table |^ we see that by extrapolating the .B-parameters from the range of masses 
accessed from our lattice to the mB d/a , without including the static HQET results, we always 
overshoot the ones that are obtained by including the static HQET values. This is especially 
pronounced for the parameter B^(m^). At this point it is not clear whether this is the real 
physical effect, or it is due to the lattice artefacts: e.g. our heavier mesons may be more 
subject to 0(a) effects, our HQET data are only perturbatively renormalized etc. Further 
research is needed to clarify this issue. 



3.2 



Physical results for (rrib) 



As in the previous section, we first convert our directly computed B^ 11 ,lovl 



4,5 



choose the MS (NDR) scheme of ref. 

(OM) 



\, according to which 



I + r 45 - 



47T 



(Ob0*)> 



RI/MOM 



55 s . We 



(38) 



16 



with the NLO matching coefficient given by |j| 

2 ( -17 + log 2 3(1 -log 2) 
-3(1 + log 2) l + log2 



^45 



(39) 



As in the previous subsection, we convert the results from table |2] at fi = 2.8(1) GeV and 
evolve them in the MS scheme to the scale /i = = 4.6 GeV. The evolution is described by 
an equation analogous to eq. Q2"T|) in which the one-loop anomalous dimension now reads [Bf 



7o 



-16 

-6 2 



whereas the NLO part (rip = 0) is 
M(/J) 



T , rMS <*s(fj) 
1 + J ' 5 47T 



/ 24379 5895 \ 
1 ~~5EUW T936" * 



J 



MS 
45 



V 



45 
T6" 



5807 
'58TJ8 



(40) 



/ 



Our 5 4 5 -parameters, in the MS scheme and at fi — mb, are given in table The matching 





q=up/down 


q= strange 




0.125 


0.122 


0.119 


0.125 


0.122 


0.119 


m P [GeV] 


1.75(9) 


2.02(10) 


2.26(12) 


1.85(8) 


2.12(9) 


2.36(10) 


Bf^(m b ) 
i?P(m 6 ) 


1.098(22) 
1.235(25) 


1.074(25) 
1.308(31) 


1.126(21) 
1.442(30) 


1.119(16) 
1.288(20) 


1.134(17) 
1.390(24) 


1.131(13) 
1.466(21) 



Table 6: B-parameters, extracted from our lattice QCD data and converted to the MS (NDR) 
scheme of ref. ^j. 



onto the corresponding operators in HQET is made through [|Kj 



n ( \ tt i (i) ®s(m P ) 
C A5 {m P ) = I + c\ 5 ' — 

47T 



45 



1/17 -11 
2 1 7 -21 



(41) 



As for the evolution of these operators in HQET, the matrix of the leading order anomalous 
dimension coefficients is 

' 7 3 



7o 



3 7 



while the NLO contribution {np — 0) reads |T0 



m = i + 



7ms 

J 45 



/ 


833 
726 + 


74vr 2 
297 


1987 
2178 


27T 2 

99 


\ 




1987 


2tt 2 


833 


74vr 2 




V 


2178 


99 


726 


+ 297 


/ 



(42) 
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We use all the above formulae to create the quantities $ 4 5 (mp, m&) ( fTT)| ) and extrapolate in 
the heavy meson mass by using eq. (|T7|). The results are 



^i,b{m Bi ,m b ) 



1.189(70) 
2.144(108) 



1.184(25) 
2.158(60) 



(43) 



After incorporating the values for _B 4 5 (mb) given in table the fit to eq. (|18|), gives 



®i${™>B d ,m b ) 



1.112(27) 
2.072(49) 



1.121(12) 
2.101(35) 



(44) 



The two ways of reaching the physically relevant results are shown in fig. |^. As in the 
previous subsection, for the final results we need to match back onto the QCD. These values 
are presented in table [7|. From tabs. ^| and |7| we conclude that, for all the £>-parameters, 



1.2 

1.0 

0.8 
2.5 

2.0 

1.5 

1.0 



i- 



0$Jm p ,m b )[\] 
•BJmJ 



I 



0$J mp ,mJ [2] 
•BJmJ 



•BL. 



0.0 



0.5 



1.0 



1.5 



1/M T 



Figure 5: The same as fig. [| but for &4^(mp, m b ) 



the extrapolations in which we do not include the static results lead to higher values. 



4 Systematic uncertainties 

The central values of this work are the ones obtained by combining the static HQET and 
QCD results, which are given in tabs. [| and [?]. 

We now need to attribute systematic errors to these results. We do not discuss the 
errors due to the use of the quenched approximation and refer to our results as quenched. 
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unconstrained 


constrained 


Light quark 


d 


s 


d 


s 


5f(m 6 ) 


1.228(72) 


1.223(26) 


1.148(28) 


1.157(13) 


5f(m 6 ) 


1.784(90) 


1.798(50) 


1.724(41) 


1.750(29) 



Table 7: Final results for the last two B-parameters defined in eq. in the MS (NDR) scheme 
of ref. j^j at \x = m&. 



The results of ref. [28 , however, are quite encouraging in that they indicate that the values 
of the parameters -81,2,3(^6) remain practically unchanged after switching from rip = to 
rip = 2. 

o Systematic uncertainties present in the static results were estimated to be in the range 
of (3 for all the .B-parameters. This error is almost entirely due to the choice of 

the renormalization point q* at which we used the boosted perturbative expressions. 
We varied 1/a < q* < it /a, and then evolved the resulting (continuum) B-parameters 
from fi = q* to fi = m^. The spread of values B(mb) with respect to the central one 
(obtained from q* = 2.6/a), has been assigned to the systematic error. 

o The QCD values are obtained after the non-perturbative renormalization in the RI/MOM 
scheme at \i — 2.8(1) GeV. We repeated the whole procedure described in the pre- 
vious section, but starting from our results obtained at /i = 1.9(1) GeV and at 
fi = 3.9(2) GeV. The final results get modified as follows: 

- AB 1 / B^ entTsl < ±1%; 



- AB 2 / B™ ntial ~ ±1 

- AB Z /Bf ntral ~ ±8 

- AB 4 /Bl cntral ~ -3%; 



- AB 5 /^ ontral ~ +11%. 

o Interpolation/extrapolation in the light quark mass is made linearly. For the average 
up/down quark mass we need to account for the possibility of the quadratic term in 
this extrapolation. As it can be seen from fig. |^, the extrapolations are smooth for all 
the bag parameters. If we include a quadratic term in the extrapolation, we obtain 
results which are fully compatible with the ones presented here. This is true for both 
QCD and static HQET 5-parameters. 

o The value a~ l {m K *) = 2.72(13) GeV, has been used throughout the paper. Another 
option would be to use the kaon decay constant, from which we obtain a' 1 (fx) = 
2.69(16) GeV. Being completely consistent with a~ 1 (mx*), this choice affects our final 
results by only +1%. 
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o Even though we use the improved action, our operators are not improved. Therefore, 
our results for the bag parameters suffer from O(a) discretization errors. The hope is 
that these errors cancel in the ratios (|5]||) from which the bag parameters are actually 
extracted. A conservative estimate on the size of these uncertainties can be obtained 
if we improve the axial current in the ratio R±, as Ao(t) — > Ao(t) + CA{P§{t + 1) — 
P 5 (t — l))/2, with the known value for the parameter ca = —0.04 ||29|| . From this 
exercise we conclude the further increase in our final results for B\ by ~ 4%. 

As for the HQET values, we checked that our values for the Bi (obtained at (3 = 6.0) 
are indistinguishable from the ones that can be extracted from the UKQCD data at 
151. 



(3 = 6.2 



That gives us more confidence that the O(a) effects in the static HQET 



data are indeed small. 

o We used the two-loop running coupling a s (fi) by taking Aq^ ' 1 = 0.25 GeV. We tried 
to vary 

A qcd 0) by 10% (which covers all the presently available lattice estimates ||30|| ), 
and see that the final results vary in the range of ±1.5%. 



We now write our results in a fully explicit form as: 



B[ d)M \m b ) = 0.87(4)(3)(0) ft) , B[ s > M \m b ) 
Bi d)m (m b ) = 0.83(3)(3)(1)(2) , ^ s)M5 (m b ) 



B. 



(d)MS 



(m fe ) = 0.90(6)(3)(7)(2) 



B 



(s)MS 



B 



(d)MS 



(m fe ) = 1.15(3)(4) ft) (3) 



B 



(s)MS 



{m b ) 
{m b ) 



B 



(d)MS 



(m b ) = 1.72(4)(5) (±58) (3) 



B^ MS (m b ) 



0.87(2)(3)(0) ft) , 
0.84(2) (3) (1) (2) , 
0.91(3)(3)(7)(2) , 
1.16(2)(4) ft) (3) , 
1.75(3)(5) ft°) (3) 



(45) 



where, besides the first statistical errors, the following sources of the systematic uncertainty 
are being written out respectively: systematics of the calculation in the static limit of HQET, 
the error in the renormalization of 5-parameters computed in QCD, combined error due to 
the variation of a -1 and of Aq^^ ' 1 (and also due to the improvement of the axial current 
in the case of Bi). After adding all systematic errors in squares we arrive at the complete 
set of results already given in table [p. 

To be able to fully reconstruct the numbers that we presented in table 0, we also need 
to provide the reader with the formulae allowing the conversion of the parameters B2(m b ) 
and B 3 (m b ) from the MS(NDR) scheme of ref. to the one of ref. f|. This is achieved 
by using the following formula 



MS [9] 



I + 



12 7T 



-11 1 

1 5 



(Os(M)> 



MS 



11 



(46) 



which we obtained after rotating the operators Qf 2 L (/-0ms °f Te ^- Hi t° the SUSY basis (|T|). 
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5 Concluding remarks 



In this paper we computed the 5-parameters for all five AB = 2 operators. The extrapo- 
lation of the results obtained directly in lattice QCD in the region of masses mp ~ 2 GeV 
to the physically interesting mass n r iB d , B , has been constrained by using the static HQET 
result. The matching QCD <-> HQET and running in each of the two theories have been 
made by the consistent use of the perturbative expressions known at NLO. The final results 
are presented in three renormalization schemes (see table |l]). 

Our results can be improved in many ways. We combined the results of the QCD lattice 
simulations performed at (3 = 6.2 with the HQET ones obtained at (3 — 6.0. Naturally, a 
good strategy would be to do the computation at the same value of (3 in both theories, to 
vary the value of (3 {i.e. of the lattice spacing) and attempt extrapolating to the continuum 
limit. All numbers are obtained in the quenched approximation (jif = 0). An investigation 
of the sea quark effects on our quenched values by repeating the analysis performed in this 
paper with tif = 2, would be very welcome. 
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Appendix: Non-perturbative calculation of the renormalization 
and subtraction constants in the (Landau)RI/MOM scheme 

In this appendix we give the numerical values for the resulting matrices of the renormal- 
ization and subtraction constants which are obtained in the (Landau) RI/MOM scheme by 
using the method of refs. |T7|, |TB|. These are computed in the following basis of operators: 

Qi = tl^ - 7sK q j l»{l ~ 7s)g j , 
Q2 = q l ln(l - 7sK q j l»{l + J 5 )q j , 
^3=^(1 + 75)^^(1-75)?^', 



(47) 



1 



i5)q ■ 



Q 4 

Qs = 2 q " a ^ 1 ~ 75 ^ q J °^ v ~ 75jg J 

which are equivalent to those appearing in eq. (|l|) (after setting b — > q): 

Qi = 1> Q 2 = -2O5 , 

Qs = 4 , Qi = O2 , 

Q 5 = 2 + 20 3 . 



(48) 



The difference between these and the results for the renormalization constants presented in 
refs. [|T3"| , |T4" is that the present renormalization and subtraction constants are not polluted 
by the Goldstone boson contributions. To eliminate those, we applied the recipe of ref. . 
Since that paper has not been released yet, we briefly explain the main steps here. 

Starting from the 4-quark Green functions computed in the Landau gauge, with all 
momenta in the external legs equal, Gi(p) = (q(p) q{p) Qi q{p) q(p)), one constructs 
the amputated ones as 



A i (p)=[l[S- 1 (p))G i (p) 



(49) 



where S 1 (p) stands for the inverse quark propagator. 
The amputated Green functions are projected onto various Dirac structures as 



Up) 



Tr 



where P~ are suitable projectors satisfying the orthogonality relation 



rfb) 



Tr 



Oij , 



(50) 



(51) 



where A.f\p) stands for the tree level amputated Green functions. The explicit ex- 
pressions for the projectors Pj can be found in ref. [I7j (eq. (37)). 
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Eq. (51) is turned into the RI/MOM renormalization condition as 



fi(p/ju) 



Sij , (52) 



where the renormalized amputated Green function, T^pf/j,), is expressed as 

fi(p//i) = r k (p)(s k j + ^{glYjZji^gl) , (53) 

up to an overall wave function renormalization, which is trivial to compute after 
imposing the vector Ward identity on the quark propagator. The condition ([52]) is 
applicable for virtualities Aq CD <C p 2 <C (n/a) 2 . 



Thus, for each of the operators from the basis (|47D, one obtains 5 equations from which 
the subtraction (A(g 2 )) and renormalization (Z^g^)) constants are computed. In 
matrix form, the final result writes 

Q(ji) = Z{n-g 2 ) [l+ A(g 2 )] Q(g 2 ) . (54) 

The structure of the Z-matrix is determined by the chirality, i.e. Qi(fi) does not 
mix with any other operator, | Q2,3(/-*) rnix with each other but not with the other 
operators. The same goes for the Q±$(ia) operators. The remaining elements of the 
5x5 matrix are filled by the subtraction constants Aij(g^). 

In practice, however, the above procedure is implemented by computing fj(p, K q ) at 
several values of the (light) quark mass (i.e. various K q ), followed by the extrapolations 
of each Aij(g%) and Zij(fi, g 2 ) to the chiral limit. This extrapolation can be dangerous 
because the operators are inserted at zero momentum (all external legs in the Green 
function have the same momentum), and the coupling to the Goldstone boson con- 
taminates the short distance behaviour (which we are interested in). In particular, 
we find that for the vertices of the structure 75 75 this coupling is indeed large. In 
addition, via projections (|50|), it may pollute the extraction of the renormalization 
and subtraction constants for the other operators. Therefore, for each projected am- 
putated four-quark Green function, one should subtract the Goldstone contribution. 
For the parity even operators (which are the ones that we consider in this paper), this 
Goldstone contribution can appear as a pole, but also as a double pole, i.e.: 

(Vi(p;K q )) =T ij (p;K q ) = Q>ij(p) + + lll^Pl 4. S i:j m 2 P . (55) 

V / j J mp rap 

Note that we also added a term 5 m 2 p , to account for the the linear dependence in the 
quark mass (m 2 p oc m,). | A judicious way to subtract the Goldstone contributions, 



5 In other words Z\\(ji) ^ 0, whereas Z\2 = Z13 = Z14 = = 0. 

6 The linear dependence in the quark mass arises after the cancellation of the quadratic quark mass term 
against the Goldstone pole contribution. More detailed discussion will be presented in ref. [pL8|. 
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and thus to reach the term a^tj?), is to consider the following combination for the fit 



with the data fL8 
m^^ijip; K qi ) - m 2 P2 T i:j ip; Kg 2 ) 7o -(p) 



a a(p) - + 5 M( m Pi + m p 2 ) ( 5e ) 

I fin III n, V / 



2 2 

m Pi ~ m P 2 '" , Pi" t P 2 

by which the pole-like contribution is automatically eliminated. This is to be per- 
formed for each value of p and accounting for all the mass combinations. [] The re- 
sulting aij{p) is thus the chiral value of the projection Tij{p) free from the Goldstone 
boson contamination. 

The procedure sketched above has been applied at three values of the renormaliza- 
tion scale: afi = 0.71, 1.03 and 1.41. In physical units, these values correspond to /i = 
1.9(1) GeV, 2.8(1) GeV and 3.9(2) GeV, respectively. The complete list of results for the 
subtraction (A^-(a)) and renormalization constants (Zij(fia)) is presented in table ||. 

As for the renormalization constants for the bilinear quark operators, which are necessary 
to compute the ratios (||) and (|[), we use the following values |3T|| : 

Zp I/MOM (l.9(l) GeV) = 0.510(9); 

^ /MOM (2.8(l) GeV) = 0.575(7); (57) 
Zf/ U0M (3.9(2) GeV) = 0.630(6). 

In addition, in the same study, we obtained Za = 0.814(4). Notice that we used the method 



of ref. |32|, to avoid the large Goldstone boson contribution to the value of ^ I / MOM g # \ s 



for the constant Za, our value is consistent with the findings of other lattice groups [^| . 



7 Besides the light masses already mentioned in the previous section, for the discussion of the renormal- 
ization constants we also had the results for n q — 0.1333 at our disposal. In this way, we could make 6 
combinations and therefore a safe extrapolation to the chiral limit. 
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Scale (fj,) 


1.9(1) GeV 


2.8(1) GeV 


3.9(2) GeV 




0.663(6) 


0.645(9) 


0.644(13) 


A 

A 12 
A13 
A14 
A15 


-0.072(3) 
-0.015(2) 
0.020(2) 
0.011(3) 


-0.069(11) 
-0.011(4) 
0.021(7) 
0.006(1) 


-0.063(10) 
-0.014(2) 
0.015(6) 
0.005(1) 


Zoo(ll) 

ZZ \r" / 


0.723(5) 


0.691(7) 


0.683(5) 


Z23 (aO 


0.315(3) 


0.257(7) 


0.202(9) 


A21 

A 24 

A25 


-0.052(2) 
-0.250(7) 
0.013(2) 


-0.055(9) 
-0.169(20) 
0.014(2) 


-0.050(9) 
-0.168(19) 
0.015(1) 




0.023(1) 


0.022(1) 


0.021(3) 


^33 (At) 


0.322(12) 


0.392(20) 


0.467(13) 


A31 
A 34 
A35 


0.018(1) 
0.351(10) 
-0.008(1) 


0.018(4) 
0.233(34) 
-0.007(3) 


0.014(3) 
0.220(31) 
-0.005(1) 


Zaa(u) 


0.414(9) 


0.473(17) 


0.534(13) 


Z45 (a*) 


-0.017(2) 


-0.015(2) 


-0.015(4) 


A41 

A42 
A43 


0.008(1) 
0.009(1) 
0.208(4) 


0.008(3) 
0.001(2) 
0.143(16) 


0.005(1) 
0.001(1) 
0.144(15) 




-0.307(4) 


-0.233(7) 


-0.176(7) 


^55 (a*) 


0.914(9) 


0.814(16) 


0.761(18) 


A51 
A 52 
A53 


0.009(1) 
0.009(1) 
0.121(1) 


0.008(2) 
0.005(3) 
0.084(6) 


0.005(1) 
0.008(1) 
0.089(8) 



Table 8: The values of the renormalization Zij(/j,) and subtraction constants Ay computed non- 
perturbatively in the (Landau ) RI/MOM scheme at three different values of the renormalization 
scale /j, at (3 = 6.2. 
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